JOURNAL OF GEOPHYSICAL RESEARCH, VOL. ???, XXXX, DOL10.1029/, 

Oxygen and hydrogen ion abundance in the 
near-Earth magnet osphere: Statistical results on the 
response to the geomagnetic and solar wind activity 
conditions 

E. A. Kronberg\ S. E. Haaland^'^ P. W. Daly\ E. E. Grigorenko^ L. M. 

4 1 5 

Kistler , M. Pranz , and I. Dandouras 



E. A. Kronberg, Max-Planck Institute for Solar System Research, Katlenburg-Lindau, Ger- 
many (kronberg@mps.mpg.de) 

S. E. Haaland, Max-Planck Institute for Solar System Research, Katlenburg-Lindau, Germany. 
Also at Department of Physics and Technology, University of Bergen, Norway. 
(stein.haaland@ift.uib.no) 

P. W. Daly, Max-Planck Institute for Solar System Research, Katlenburg-Lindau, Germany 
(daly@mps.mpg.de) 

E. E. Grigorenko, Space Research Institute, Russian Academy of Sciences, Moscow, Russia 

(elenagrigorenko2003@yahoo . com) 
L. M. Kistler, University of New Hampshire, Durham, NH, USA (Lynn.Kistler@unh.edu) 
M. Franz, Max-Planck Institute for Solar System Research, Katlenburg-Lindau, Germany 

(fraenz@mps.mpg.de) 

I. Dandouras, IRAP, University of Toulouse, UPS-OMP, CNRS, Toulouse, Prance (Ian- 
nis. Dandouras@irap.omp.eu) 



DRAFT 



October 17, 2012, 12:22ain 



DRAFT 



X - 2 KRONBERG ET AL.: OXYGEN ABUNDANCE IN THE MAGNETOSPHERE 

5 Abstract. The composition of ions plays a crucial role for the fundamen- 

6 tal plasma properties in the terrestrial magnetosphere. We investigate the 

7 oxygen-to-hydrogen ratio in the near- Earth magnetosphere from -10 Re < 

8 Xqse < 10 Re- The results are based on seven years of ion flux measure- 

9 ments in the energy range ~10 keV to ~955 keV from the RAPID and CIS 

10 instruments on board the Cluster sateUites. We find that (1) hydrogen ions 

11 at ~10 keV show only a slight correlation with the geomagnetic conditions 

12 and interplanetary magnetic field changes. They are best correlated with the 

13 solar wind dynamic pressure and density, which is an expected effect of the 

14 magnetospheric compression; (2) ~10 keV 0+ ion intensities are more strongly 
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15 affected during disturbed phase of a geomagnetic storm or substorm than 

16 >274 keV 0+ ion intensities, relative to the corresponding hydrogen inten- 

17 sities; (3) In contrast to ~10 keV ions, the >274 keV O"*" ions show the strongest 

18 acceleration during growth phase and not during the expansion phase itself. 

19 This suggests a connection between the energy input to the magnetosphere 

20 and the effective energization of energetic ions during growth phase; (4) The 

21 ratio between quiet and disturbed times for the intensities of ion ionospheric 

22 outflow is similar to those observed in the near-Earth magnetosphere at >274 

23 keV. Therefore, the increase of the energetic ion intensity during disturbed 

24 time is more likely due to the intensification than to the more effective ac- 

25 celeration of the ionospheric source. In conclusion, the energization process 

26 in the near-Earth magnetosphere is mass dependent and it is more effective 

27 for the heavier ions. 
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1. Introduction 

28 There are two major sources for the ion population in Earth's magnetosphere: the solar 

29 wind and the ionosphere. 

30 Hydrogen, helium and oxygen are the most abundant elements in the solar wind. Ac- 

31 cording to measurements by SWICS on Ulysses, the oxygen abundance (H/0 ratio) varies 

32 significantly between fast and slow solar winds, with mean values of 1500 and 2300, re- 

33 spectively [e.g., von Steiger et al, 2010]. A direct entry of solar wind plasma into the 

34 magnetosphere is believed to occur mainly when the interplanetary magnetic field (IMF) 

35 is oriented southward. This leads to reconnection first at the dayside magnetopause and 

36 then in the magnetotail. Also, northward IMF results in reconnection at southern and 

37 northern high latitudes, facilitating the entry of solar wind plasma. 

38 The main source of oxygen in the magnetosphere is believed to be the terrestrial iono- 

39 sphere [Yau et al, 1984; Lockwood et al, 1985; Gloeckler et al., 1985; Chappell et al., 

40 1987; Yau and Andre, 1997; Chappell et al., 2000; Cully et al., 2003b; Huddleston et al., 

41 2005; Moore and Horwitz, 2007; Kitamura et al, 2010], where oxygen ions are accelerated 

42 upward by electric fields parallel to the background magnetic field and pressure gradients. 

43 The ionospheric ions are either accelerated in the cusp/cleft region and then convected 

44 across the polar cap into the lobe [Lockwood et al, 1985; Kistler et al., 2010], or they 

45 come from the nightside aurora, which provides a fast feeding of O"*" ions in the near- 

46 Earth plasma sheet during the substorm expansion phase [Yau et al., 1985; Daglis and 

47 Axford, 1996]. The kinetic energy of upfiowing ions is typically a few tens of eV [Yau 

48 and Andre, 1997], whereas typical plasma sheet energies are much higher [e.g., Lennarts- 
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49 son and Shelley, 1986]. Centrifugal forces accelerate particles as they travel tailward 

50 along field lines [e.g., Cladis, 1986; Nilsson et al., 2008, 2010], but the most significant 

51 energization will take place after the particles enter the plasma sheet. 

52 The accelerated plasma sheet particles are convected Earthward and populate the ring 

53 current which influences the transport in the radiation belts. The radiation belts are a key 

54 region for space weather studies, as their dynamics can induce damage to satellites and 

55 other space weather effects. Forecasting models based on radiation belt dynamics include 

56 the dynamics of the ring current. Additional ionospheric particles affect the background 

57 magnetic field, densities, temperature and convection electric field, and therefore lead to, 

58 e.g., changes in the ring current [Welling et ai, 2011]. It is challenging to include the 

59 heavy ion contribution in numerical studies of the ring current [e.g., Glocer et al, 2009; 

60 Welling et ai, 2011] or radiation belt models. Knowledge about the ion composition and 

61 its dependence on solar wind and magnetospheric conditions is important for defining the 

62 boundary conditions in such simulations. 

63 A number of earlier studies have demonstrated that the abundance of oxygen in the ring 

64 current region can change significantly with geomagnetic activity. Sharp et al. [1981] have 

65 reported measurements of oxygen in the energy range 0-17 keV/e at distances between 

66 11 and 23 Re downtail by the ISEE 1 spacecraft. Gloeckler et al. [1985] analyzed for the 

67 first time the ion composition measurements of the ring current during the storm time 

68 using AMPTE/CCE spacecraft data in the energy range ~1 to 300 keV/e at L values 

69 between ~3 and 6.5. They suggested that large changes in the O"*" number and energy 

70 densities implied an injection of energetic ionospheric ions. Hamilton et al. [1988], using 

71 AMPTE particle measurements, found an event in which the O"*" dominates H"*" in the ring 
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72 current during a large magnetic storm. A later study by Daglis et al. [2000] established a 

73 connection between 0+ (50-426 keV, measured by CRRES) enhancements during storms 

74 and substorm activities in the inner magnetosphere at L~7. Another, more extensive 

75 study by Fu et al. [2002], using CRRES observations, has shown that not all injection 

76 events have the same ion composition. Their statistical results demonstrated that about 

77 73% of all the events were rich in oxygen ions and the rest have very low admixtures 

78 of oxygen. Events without oxygen ions were found to be mainly correlated with weak 

79 geomagnetic activity conditions. 

80 The abundance of oxygen in the magnetotail's plasma sheet can vary significantly with 

81 geomagnetic activity. Moebius et al. [1987] showed that the spectra harden in the mag- 

82 netotail's plasma sheet during substorm, with a larger energy increase for 0+ than for 

83 H"*", i.e., the 0"'"/H+ ratio also increases with energy. For this case study they used ob- 

84 servations from AMPTE/IRM which resolved the ion composition of the suprathermal 

85 population (10-230 keV/e). Kistler et al. [1990] used AMPTE/IRM observations to com- 

86 pare spectral changes at 15-19 Re (tail) and 7-9 Re (inner magnetosphere, night side). 

87 They found that the spectral changes were very similar at both locations. Kistler et al. 

88 [1994] investigated how the flux of the energetic particles observed by the AMPTE/IRM 

89 decreases after the initial particle "injection". These observations show that the decay 

90 rate is organized by total energy, energy per charge and velocity. In a statistical study, 

91 Nose et al. [2009] used 16 years of Geotail data in order to look at the suprathermal ion 

92 composition in the plasma sheet in the range —IOORe < Xqsm < —8Re- They investi- 

93 gated the correlation between the 0+/H+ ratio and solar activity using the F10.7 index 

94 and with the geomagnetic activity using the Kp index. They found that physical pro- 
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95 cesses in the plasma sheet are expected to be much different during solar minimum and 

96 solar maximum, because the Alfven velocity changes significantly between these two solar 

97 activity regimes. More recently, ion composition data from the CIS/CODIF instrument 

98 on the Cluster spacecraft were used to study 0+ and H+ bulk content (0-40 keV/e) in 

99 the plasma sheet within 15-19 Re as a function of solar activity and geomagnetic activity 

100 [Mouikis et al., 2010]. They found that the 0+ density were strongly dependent on these 

101 parameters. 

102 In this paper, we present a comprehensive study of the O"*" and H+ abundance in the 

103 plasma sheet of the near-Earth magnetosphere (from -10 Re < Xqse < 10 Re)- The 

104 data are obtained from seven years of particle measurements by the Cluster satellites. 

105 We compare the characteristic behavior of the ions at ~10 keV and at 274-955 keV to 

106 different solar wind and geospace activity factors. 

107 The novelty of this study is a combination of (a) the extensive spatial region coverage of 

108 the near-Earth magnetosphere, (many previous studies either focused on the magnetotail 

109 region or were mainly case studies in the near-Earth region only); (b) energy range, as 
no no extensive statistical study has been done so far with energies up to ~955 keV; (c) 

111 a comparative analysis of ion intensities at different energies, namely at ~10 keV and 

112 from 274 keV to ~955 keV; (d) the advantage of having more accurate parameters of the 

113 solar-terrestrial coupling which became available during the last 15 years. 

114 We can for the first time investigate the details of how the near-Earth hydrogen and 

115 oxygen populations depend on the geomagnetic and solar wind activity. 

116 The paper is organized as follows: In Section 2, we give a brief overview of the Cluster 

117 spacecraft, the RAPID and CIS ions spectrometers and their data products. In Section 3, 
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118 we explain how we process the data and derive the 0"'"/H+ ratio. Section 4 presents the 

119 observations and shows how the 0+ and H+ abundance varies with different parameters 

120 controlling the dynamical state of the magnetosphere such as solar wind conditions and 

121 also geomagnetic disturbance level. Section 5 discusses the results and related physical 

122 processes. Section 6 summarizes the observations. 

2. Data and instrumentation 

123 The results presented in this study are primarily based on in-situ measurements from 

124 the Cluster quartet of spacecraft. The Cluster mission comprises four identical spacecraft 

125 flying in a tetrahedron-like formation. More information about the Cluster mission and 

126 instrumentation is given in Escoubet et al. [1997]. In this paper we have used data from 

127 the SC4 spacecraft, since this gives the best data return for our purpose. 

128 For the period covered by this study Cluster was in a nearly 90° inclination elliptical 

129 orbit with apogee around 19.6 Re and perigee around 4 Re. Due to spacecraft orbits in 

130 the near-Earth region our observations mainly cover the dawn and dusk flanks and the 

131 dayside plasmasheet, mainly the north part, see Figure 1. This Figure shows maps of the 

132 energetic O"*" observations, projected into XYqse, YZose dayside and YZqse nightside 

133 planes, respectively. This coverage allows us to study the transport of ions from the tail 

134 region around the Earth. 

135 In addition to the in-situ Cluster observations, we also use solar wind magnetic field and 

136 plasma data and geomagnetic disturbance indices to check dependencies and correlations. 

137 In this study, we use the AE and Dst indices and not the Kp index as e.g. Nose et al. 

138 [2009] or Mouikis et al. [2010] because they are more precise (the Kp index has a 3 hours 

139 cadence, while the AE and Dst indices have higher time resolution) and have a more direct 
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140 physical meaning. The Dst index primarily reflects slow (hours) magnetic perturbations 

141 caused mainly by ring current enhancements, whereas the AE index reflects short time 

142 (minutes) perturbations at high latitudes caused by variations in fleld-aligned current due 

143 to plasma sheet processes such as bursty bulk flow activity. 

2.1. Cluster peirticle measurements 

144 We use ion measurements from two different sensors onboard Cluster. Energies above 

145 30 keV are covered by RAPID: the Research with Adaptive Particle Imaging Detector 

146 [see Wilken et al., 2001], whereas the lower part of the energy range (energies below 40 

147 keV) is covered by CIS: the Cluster Ion Spectroscopy instrument [see Reme et al, 2001]. 

148 2.1.1. Energetic ion fluxes from RAPID 

149 The RAPID Imaging Ion Mass Spectrometer (IIMS) consists of three identical detector 

150 heads and use a combination of time-of-flight and energy measurements to classify and 

151 bin the incident particles. The three sensors cover 180° in polar angle, and by utilizing the 

152 spacecraft spin, RAPID/IIMS provides the distribution of ions with complete coverage 

153 of the unit sphere in phase space. Due to a degradation of the central heads on the ion 

154 sensor on all spacecraft during the early phase of the mission, only a limited part of the 

155 full 3D distribution is available during later stages, however. Since our purpose is to 

156 investigate the statistical abundance of two distinct ion species in a particular plasma 

157 region, we are not dependent on any 3D abilities, but integrate over all directions to 

158 get an omnidirectional flux. Except for extremely anisotropic distributions, the lack of 

159 sensitivity in the central sensor does not play any major role [Kronberg et al, 2010]. 

160 RAPID/IIMS allows for discrimination of various ion species by utilizing time-of-flight 

161 and energy information of each incident particle. The resolution of the instrument makes 
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162 it possible to group ions into three categories; protons (H"*" atomic mass unit 1), helium 

163 (He, atomic mass units 2-4) and a common group of heavy ions (atomic mass units 14- 

164 16). This latter group is dominated by oxygen, but the mass resolution does not allow 

165 for distinction between oxygen, carbon and nitrogen. This group is therefore commonly 

166 refereed to as the CNO group. The energy ranges are 28 keV - 4 MeV for protons, 137 

167 keV - 3.8 MeV for hehum and 274 keV - 4 MeV for the CNO group. No information 

168 about charge state is possible. For simplicity, we will simply refer to the ratio 0'''/H"'" in 

169 case of RAPID data, although, we mean the CNO mass group and hydrogen, respectively. 

170 The time resolution of the RAPID omnidirectional data is 1 spin (~4 seconds), but 

171 since the count rates are often very low, we have used 1 minute averages in this study. 

172 2.1.2. Moments and pcirticle measurements from CIS 

173 CIS consists of two sensors. In this paper we use data from the Composition and 

174 Distribution Function (CODIF) instrument to calculate an omnidirectional flux as done 

175 for RAPID. The energy range for CIS is 0.03-40 keV per charge. In addition to the 

176 omnidirectional proton and O"*" fluxes, we also utilize plasma moments from the CODIF 

177 sensor. 

178 As for RAPID we use 1 minute averages and put them on the same timeline as the rest 

179 of the data set. Plasma moments are used to identify the plasma region and to fllter the 

180 data. The CIS instrument has several sensitivity modes, automatically switched in-flight 

181 in order to achieve optimal sensitivity and resolution for a particular plasma region. To 

182 avoid errors introduced by mode shifting and non-optimal modes we only use data from 

183 the so-called magnetospheric mode 13 [see details in Reme et al, 2001; Dandouras et al, 

184 2006]. 
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185 0+ measurements by CIS can be slightly contaminated by protons (the major species). 

186 The contamination of H+ into O"*" is about 0.1-0.3% (depending on energy and on a 

187 particular mode) for low count rates (e.g. in the magnetosphere) for ~10 keV ions. 

188 In order to compensate for this contamination of oxygen by protons during active time 

189 periods, we calculate the oxygen intensities as follows: j*^^ — — j^^ ■ 0.001, where 

190 and j^^ are the measured oxygen and proton intensities, respectively. 

3. Methodology 

3.1. Calculation of the 0+/H+ ratio 

191 Our purpose is to compare 0'''/H"'", namely, the ratio between intensities which are 

192 proportional to the number fluxes of 0+ and H+ [1/cm^ sr s keV] in two energy ranges. 

193 This ratio is equal to the ratio between integral fluxes if the energy range is identical for 

194 the two species. To be able to identify any spectral changes, we use two different energy 

195 ranges, ~10 keV (8.1 to 10.3 keV) and the energy band 274-955 keV. These energy bands 

196 are somewhat arbitrarily chosen, and are a compromise between the energy binning of the 

197 two instruments and sufficient statistics. We take a wide energy band for the energetic 

198 ions since the measurements of ions at these energies have relatively low count rates and 

199 wide energy channel divisions. As we do not consider details of any particular acceleration 

200 process but only the response of the particle population to the different geomagnetic and 

201 solar wind activity, the use of the wide energy range of the energetic ions is justifled in 

202 this study. 

203 For the CIS instrument, the ratio between the number fluxes of O"*" and is based on 

204 the same energy range; since the energy channels are similar for protons and oxygen. For 

205 RAPID, however, the energy binning for 0+ and H+ have different thresholds, see energy 
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06 ranges listed in Table 1. The details on how the RAPID energy channels were rebinned for 

07 the most effective intensities comparison are described in the Appendix. These RAPID 

08 intensities (J^^ and J'^^) are used in the current study. 

3.2. Construction of plots 

09 Figures 2 to 6 show median intensities for protons and oxygens in two different en- 

10 ergy bands versus the respective geomagnetic or solar wind parameter. The plots are 

11 constructed from data accumulated in the plasma sheet region, defined by plasma beta 

12 value in the range 0.1-10 [Baumjohann et al, 1989]. We have excluded observations 

13 where the magnetic field component was negative in order to avoid the cusp re- 

14 gion. Furthermore, we limit our study to the near-Earth magnetosphere, i.e. distances 

15 — 10Re< Xgse < 10 Re- We also exclude data from radial distances < 6 Re, as the 

16 sensors can be affected by pileup effects from both penetrating particles (e.g. during solar 

17 proton events) and high count rates. 

18 To minimize the effect of skewed distribution, we use median rather than mean values 

19 of the distribution functions. For the same reason, statistical spreads are given as median 

20 absolute deviations (MAD), [see e.g., Venables and Ripley, 1999]. Median value and MAD 

21 are less sensitive to non-normal distribution functions than the mean value and standard 

22 deviation. In order to determine whether differences between values are significant, con- 

23 fidence interval (CI) error bars are given in Figures 2 to 6: CI — SE ■ where tn-i is 

24 the t-distribution with n degrees of freedom, SE = a j \/n is standard error, where cr, the 

25 standard deviation, was calculated in this case as: cr = 1.4826 • MAD [Huber, 1981]. 

26 The 0"'"/H+ ratio is defined as the ratio between median values of the O"*" and H+ 

27 intensities. Correspondingly, errors for the 0+/H+ ratio are calculated as: CIo+/h+ — 
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228 {CIo+ /0~^Y + {CIh+ / H+y -O^ /Yi'^ . Each value shown in the Figures contains at least 

229 100 measurements. 

230 The histogram ranges for Dst, AE, solar wind density and solar wind pressure are de- 

231 fined as follows. The distributions of these parameters are skewed. In order to avoid 

232 outliers, which would take attention away from the more typical values we (1) compute 

233 quartiles of the distribution; (2) define interquartile range (IQR) between first quartile, 

234 Ql, and third quartile, IQR=Q3-Q1; (3) use the 1.5-IQR rule for outliers, which approxi- 

235 mately corresponds to 3cr. This will lead to the histogram ranges between Q1-1.5-IQR and 

236 Q3+1.5-IQR. The behavior of the oxygen and proton intensities are rather complicated 

237 at values which are larger than ±3cr. As an exception we extended the range for the Dst 

238 from -70 nT to -100 nT, in order to cover the whole range for moderated storms. 

4. Results 

239 Below, we present and discuss how the 0+, H+ intensities and 0+/H+ ratio depend on 

240 respective parameters of the solar-terrestrial coupling in more detail. 

4.1. Response to magnetic activity 

241 4.1.1. Geomagnetic storms 

242 The interaction between a solar wind shock wave and/or an interplanetary magnetic 

243 cloud with the terrestrial magnetic field can lead to the phenomena called magnetic storms 

244 during which the horizontal component of the Earth's magnetic field decreases dramati- 

245 cally. This decrease is reflected by the Dst index. 

246 The response of the ion intensities to the changes in the Dst index is shown in Figure 

247 2 and summarized in Table 2. The intensity of protons at ~10 keV and at >274 keV 
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248 increases approximately in the same manner with decrease of Dst index down to -100 

249 nT (the maximum intensity is ~2 times higher than the minimum intensity). Oxygen 

250 intensities at both energies show more dramatic changes with Dst index. Intensities are 

251 up to 30 times higher between the minimum to the maximum. However, during weak 

252 (Dst below -30 nT) and moderate (Dst below -50 nT) magnetic storms, the ratio O^/H"*" 

253 is affected more strongly at ~10 keV. Since the ~10 keV 0+ intensity grows significantly 

254 (5 times) at Dst from -30 to -100 nT, compared to no change of 0+/H+ ratio at >274 

255 keV. At the higher energies, the intensities of both species grow together during weak 

256 and moderate storms, therefore, the 0+/H"'" ratio does not clearly change. This agrees 

257 with Ono et al. [2009] who found that during dipolarization associated with substorms the 

258 energy density of energetic oxygen does not always increase. However, we see a significant 

259 growth of the >274 keV oxygen, and of 0"'"/H"'" ratio during the growth phase, from 15 

260 to -30 nT. Therefore, a significant amount of >274 keV oxygen is effectively accelerated 

261 during a storm growth phase. 

262 At positive Dst>15 nT the intensity of >274 keV protons is clearly higher compared to 

263 that at Dst~ 0. This is expected as positive Dst values are associated with a compression 

264 of the magnetopause which can lead to a drift shell displacement. An enhanced intensity 

265 of ~10 keV H"*" is observed for positive values of Dst. The effect is not observed in the 

266 0+ intensity, though. 

Based on the observations, we can define linear functional dependence between the 
O^/H"*" ratio and disturbance level (from -60 to 30 nT): 

0+/H+ = -9.7 • 10"^ X Dst[nT] + 3.5 • 10"^. 

267 The dependency is constructed using a linear least squares fit. 
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268 We find a weak correlation between proton intensity (both at ~10 keV and >274 keV) 

269 and Dst. This is consistent with results reported by e.g. Nose et al. [2009]; Mouikis et al. 

270 [2010]. Also consistent with our results, Mouikis et al. [2010] found significantly stronger 

271 correlation of oxygen intensities and geomagnetic activity in the magnetotail region. 

272 4.1.2. Response to magnetotail processes 

273 Processes in the magnetotail due to e.g. magnetospheric substorms and bursty bulk 

274 flows causing enhanced field-aligned currents leading to auroral electrojet activity reflected 

275 in the AE index. 

Figure 3 shows the response in H+ and 0+ intensities to these processes. The increase 
in the AE index is associated with an exponential growth by a factor 4 (from the minimum 
to the maximum) of the ~10 keV oxygen intensities. The ~10 keV protons only show a 
slight increase with AE index, (maximum value ~1.5 time higher than the lowest value). 
As above, we can now establish an empirical functional relation between 0"^/H"^ ratio 
and AE index for ~10 keV ions 

0+/H+ = 8 • 10"^ X AE[nT] + 3.5 • 10~^ 

276 The increase of the O'^/R'^ ratio at ~10 keV energies means that the oxygen is accelerated 

277 up to these energies more effectively than protons, although sources of their acceleration 

278 have to be different. The study by Ono et al. [2009] shows similar results: acceleration 

279 of 0+ at energies of 9 to 36 keV/e during substorms (by electric fields induced by mag- 

280 netic field fluctuations whose frequencies are close to the gyrofrequencies of ions) is more 

281 effective than of H+. Event studies by Zong et al. [2008, 2009] have shown that bursty 

282 bulk flows can be responsible for supply of the ionospheric oxygen into the near-Earth 

283 magnetosphere. 
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284 Proton intensities at >274 keV show an exponential increase with AE. Intensities are 

285 ~8 times higher for high AE values. A similar dependence between AE and >274 keV 

286 oxygen intensities are also observed, although the slope has a break at 150-250 nT (strong 

287 growth before and moderate after). The highest intensities of >274 keV oxygen, at high 

288 AE values are about 33 times higher than the lowest intensities. 

289 The effect of substorms on the abundance of the 0+ ions relative to protons is larger 

290 at ~10 keV, since the O"*" intensity increases more strongly than H+. At higher energies, 

291 the intensities of the two species seem to increase approximately in the same way starting 

292 from AE~200 nT. This again agrees with results obtained by Ono et al. [2009] who did 

293 not find a clear correlation between substorm-associated dipolarizations and effective 0+ 

294 acceleration for ions in the energy range 56-212 keV/e. However, our results show a 

295 significant increase in intensities of >274 keV oxygen from quiet (AE~0-100 nT) to more 

296 disturbed times (AE~150-250 nT). Therefore, a significant amount of >274 keV oxygen 

297 is effectively accelerated during weak/moderate substorms or during a substorm growth 

298 phase. 

4.2. Response to solar wind conditions 

299 4.2.1. Response to IMF direction 

300 The response of the ion intensities to the changes in IMF direction is shown in Figure 

301 4. The IMF direction is represented by the clock angle (CA) defined as the angle between 

302 Zqse and the projection of IMF into the YZqse plane. 

303 When IMF Bz is negative (i.e. southward IMF), the couphng to the Earth's geomagnetic 

304 field is strongest. Southward interplanetary magnetic conditions are typically associated 

305 with geomagnetic storm activity. 
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06 In order to take into account that the changes in particle population do not immediately 

07 respond to changes in IMF direction, we only use those measurements during periods 

08 where IMF was reasonably stable for at least 30 minutes. 

09 The intensity of ~10 keV protons increases slightly at the IMF clock angle, -115° in 

10 ~1.5 times from the minimum to the maximum. For the oxygen intensities the increases 

11 are observed at -115° and 125°, in ~2.5 times from the minimum to the maximum. 

12 High energy particles (energy >274 keV) are much more sensitive to IMF clock angle 

13 variations, especially oxygen ions. The oxygen and proton intensities at >274 keV are 

14 a factor of 10 and 3 higher for southward IMF, respectively. We find a strong positive 

15 correlation between the intensities of >274 keV O"*", H+, the O^/H"*^ ratio and the IMF 

16 direction. The energetic ion intensities seem to be correlated with the "openness" of the 

17 magnetosphere or change of the clock angle towards the southward direction. However, 

18 during the most southward IMF (and dawnward) the acceleration ceases to be effective, 

19 the O^/H"*" ratio does not show significant changes. 

20 Since RAPID does not discriminate between charge states, oxygen of solar wind origin 

21 (with higher charge states) might explain the strong correlation with IMF clock angle. 

22 To check this, we compared our intensities with corresponding intensities of high charge 

23 state oxygen, with clear solar wind origin, measured by Kremser et al. [1987] in the inner 

24 magnetosphere using the AMPTE/CHEM data set. Their results indicate that intensities 

25 of high charge state oxygen were about one order of magnitude lower than our RAPID 

26 results. Hence, we think that our observations are primarily of ionospheric origin. The 

27 effect of IMF orientation is significantly larger for the intensity ratio of the 0+ ions to 

28 protons at >274 keV than at ~10 keV. These all suggest a connection between the energy 
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329 input to the magnetosphere related to the change of clock angle and the energization of 

330 0+ to high energies (>274 keV). 

331 4.2.2. Solar wind density and pressure 

332 The response of the ion intensities to changes in solar wind density and pressure is 

333 shown in Figures 5 and 6, respectively Intensities of ~10 keV oxygen and protons are 

334 almost exponentially increasing with solar wind density, (by a factor of ~4 for H+ and a 

335 factor of ~3 for O"^. For energies >274 keV, we do not see any clear dependence on the 

336 solar wind density 

337 The ~10 keV protons and oxygen show almost exponential increase (by a factor ~ 4) 

338 with the solar wind dynamic pressure. The proton intensities at >274 keV increase (by a 

339 factor of ~ 5) with the solar wind pressure. There is no clear general dependency of the 

340 >274 keV oxygen intensities on the solar wind pressure. However, for >274 keV O"*", we 

341 observe a ~40 fold increase between ~0.5 nPa to 1.75 nPa. At these solar wind pressure 

342 values, the energetic oxygen is accelerated more effectively than energetic protons. This is 

343 also supported by the fact that in ~10 keV oxygen such a tendency is not observed. From 

344 ~1.75 nPa to 3nPa the intensity of the >274 keV oxygen and protons becomes almost 

345 equal. 

346 Statistically, the solar wind density strongly correlates with the plasma sheet density 

347 [Borovsky et ai, 1998]. This is seen for the ~10 keV protons. For ~10 keV protons we 

348 find a strong correlation between intensity and the solar wind density and pressure. 

349 The increase of ~10 keV O"*" intensities is well correlated with the increase of the solar 

350 wind dynamic pressure. This is in agreement with earlier studies, e.g. by Cully et al. 

351 [2003a] ; Nose et al. [2003] who found that the ionospheric outflow is strongly affected by 
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352 the solar wind dynamic pressure. Also in a study by Echer et al. [2008], it was shown that 

353 a shock compression can lead to enhanced inflow of the oxygen ions into the near-Earth 

354 magnetosphere. However, this can be also compression effects, see Section 5.5. 

355 Intensities of ~10 keV oxygen ions show clear dependencies on solar wind dynamic 

356 pressure, as this energy is closer to source energies. This is not the same for energetic ions 

357 which are not closely related with outflow energies and require stronger acceleration. The 

358 observed prominent increase in the oxygen intensities for solar wind dynamic pressure up 

359 to 1.5 nPa is more hkely related to acceleration mechanisms due to the thinning of the 

360 plasma sheet, see Section 5.4. The largest increase in the energetic particle intensities is 

361 seen during the transition period between quiet and disturbed times. 

5. Discussion 

5.1. Compcirison to other observations 

362 Our energy density results listed in Table 3 (see calculations in Appendix, Equation ) 

363 are comparable to the observations in the near-Earth plasma sheet made by e.g. Gloeckler 

364 et al. [1985], Hamilton et al. [1988] and Nose et al. [2001]. We do not observe the dramatic 

365 differences between quiet and disturbed times reported in these earlier studies, though. A 

366 possible reason is that we have used median values rather than mean values, and thus avoid 

367 effects of outliers and strongly scattered distributions. Also, the definitions of quiet and 

368 disturbed times, are based on AE index and Dst indices and differs from earlier studies. 

369 One could notice that the difference between quiet and disturbed times for the ~10 keV 

370 ions is stronger when using the Dst index as definition. Taking even more negative values 

371 of Dst will lead to the higher 0+/H+ ratios. Therefore, the dynamics of the ~10 keV 

372 particles in the near-Earth magnetosphere is more strongly correlated with the Dst index. 
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5.2. Comparison to ionospheric outflow observations 

373 Ion outflow rates (at 0.01-17 keV) can vary by a factor of 30 between quiet and very 

374 disturbed periods Yau et al. [1985]. Our result shows that the >274 keV oxygen intensity 

375 is about 33 times higher during very disturbed periods than during quiet time (here AE 

376 index is used). For comparison, the ~10 keV O"*" intensity increases a factor 4 between 

377 quiet and very disturbed times. Acceleration of ~10 keV oxygen seems to be more effective 

378 during disturbed times, while for >274 keV oxygen acceleration primary occurs during 

379 transitions from quiet to disturbed periods. Such slope break between quiet and disturbed 

380 states is not observed in the ionospheric 0+ outflow [Yau et al., 1985; Cully et al., 2003a]. 

381 For protons, the intensities of ~10 keV ions are only slightly different during disturbed 

382 periods than during quiet times. The intensity of protons with energies >274 keV increases 

383 by a factor of about 8 between quiet and disturbed times. For comparison, for auroral 

384 and polar protons outflow an increase of the order of 5 was reported by Yau et al. [1985]. 

385 The ratio between quiet and disturbed times for the ion ionospheric outflow is higher 

386 than those observed in the near-Earth magnetosphere at ~10 keV but similar to those 

387 observed at >274 keV. These appears to be a connection between the energetic ion inten- 

388 sities and ion outflow that changes with geomagnetic activity. The fact that the ratios 

389 between quiet and disturbed times are higher for the more energetic and heavier particles 

390 means that the ions are more effectively accelerated up to these energies. This is consis- 

391 tent with explanation by Moebius et al. [1987] and many other studies thereafter, that 

392 the acceleration process in the near-Earth magnetosphere is mass dependent, as it is more 

393 effective for the heavier particles. 
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5.3. Response of protons to geomagnetic activity and solar wind conditions 

394 As shown in Section 4.2.1, the intensity of ~10 keV protons does not seem to be directly 

395 affected by dayside reconnection. Response to IMF direction would be expected, though, 

396 as it is likely to be associated with dayside reconnection and then consequent transfer of 

397 the solar wind ions into the magnetosphere and subsequent acceleration. By the same 

398 token, magnetic storms and substorm activity are typically also related to the dayside 

399 reconnection. In addition storms and substorms are typically associated with enhanced 

400 outflow. As a consequence, more ions would be energized in the magnetosphere. Still, the 

401 correlation between the ~10 keV proton intensities on one hand, and geomagnetic activity, 

402 solar wind dynamic pressure and clock angle on the other hand is less pronounced than 

403 for the ~10 keV oxygen, see Sections 4.1.1, 4.1.2 and 4.2.1. 

404 These observations can be explained as follows. Heavier ions (0+) require additional 

405 forces such as electric fields, wave activity or stronger gradients in pressure and tempera- 

406 ture to escape from the ionosphere. These processes are primarily driven by enhanced day- 

407 side reconnection and the subsequent processes reflected by increased disturbance levels. 

408 Protons, on the other hand, with their much lower escape energies, are not so dependent 

409 on these external forces. Solar illumination alone suffices to extract a significant number 

410 of light ions from the ionosphere. Therefore, the strong dependence on the disturbance 

411 parameters of the ~10 keV protons is not observed. 

5.4. Changes during growth phase 

412 Why would we expect an increase in energetic oxygen intensities during the transition 

413 period between quiet and disturbed times (substorm growth phase)? The increased solar 

414 wind pressure (or loading of the energy during Dungey cycle) can lead to the plasma sheet 
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415 thinning [e.g., Sauvaud et ai, 1996] and therefore to a more effective acceleration of the 

416 heavy ions [e.g., Ganguli et al, 1995]. There are different regimes of particle motion in the 

417 plasma sheet derived by Biichner and Zeleny [1986]. According to their work the particles 



418 are accelerated most effectively when parameter k — ^J^^f^ ~ 1, where Vcurv is the radius 



419 of magnetic field curvature and Vg is ion gyroradius. Under these conditions the particle 

420 trajectories become stochastic and particles are effectively accelerated by quasi-stationary 

421 dawn-dusk electric fields. A decrease in plasma sheet thickness, e.g. in response to an 

422 increase of the lobe magnetic field, will lead to larger region where stochastic acceleration 

423 of energetic ions can take place, as it will decrease the parameter k which is quite large 

424 in the plasma sheet for solar wind protons [Ashour-Abdalla et al, 2009]. 

5.5. Changes with solar wind pressure 

425 Increased solar wind dynamic pressure leads to a compression of the magnetosphere and 

426 the same total quantity of particles have to be distributed over a smaller volume. Con- 

427 sequently, the particle density /intensity will be higher and this will lead to the observed 

428 correlation. The observed differences in behavior between ~10 keV (where a correlation 

429 is evident) and > 274 keV ions (no clear correlation with dynamic solar wind pressure) 

430 can be explained as follows: the plasma sheet may become thinner as result of very strong 

431 solar wind pressure, and therefore, the gyroradius versus curvature radius effect discussed 

432 above may be relevant also here. 

433 On other hand, the correlation with the solar wind dynamic pressure can be related to 

434 the Earth's ionosphere which may provide heavier ions to the magnetotail during enhanced 

435 solar wind dynamic pressure, see Cully et al. [2003a]; Nose et al. [2003]; Echer et al. [2008]. 
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6. Summary 

436 Based on seven years of ion composition data from Cluster observations at radial dis- 

437 tances -10 Re < Xqse < 10 Re, we find the following: 

438 (1) H"*" intensities at ~10 keV show only a shght correlation with geomagnetic conditions 

439 (Dst and AE indices) and interplanetary magnetic field orientation. Solar wind dynamic 

440 pressure and density (i.e. effects of magnetospheric compression) seem to have a larger 

441 effect on the ~10 keV H"*" intensities. 

442 (2) 0+ ion intensities at ~10 keV are more affected by geomagnetic storms and sub- 

443 storms than >274 keV ion intensities, than their corresponding hydrogen counterparts. 

444 The >274 keV O"*" energization/acceleration seem to be strongest during the transition 

445 period from quiet to disturbed times, i.e., during growth phases rather than during the 

446 disturbed phases itself. The ~10 keV ions do not reveal such a dependence. 

447 (3) We find a strong positive correlation between the flux of >274 keV O"*", H+, the 

448 0"'"/H"'" ratio and the IMF direction. This demonstrates a connection between energy 

449 input to the magnetosphere and effective energization of energetic ions. 

450 (4) The intensity ratio between quiet and disturbed times for the ionospheric ion outflow 

451 is similar to those observed in the near-Earth magnetosphere at >274 keV. Therefore, the 

452 observed increase of the energetic ion intensity during disturbed time is not only due to 

453 a more effective acceleration but also due to enhanced ion outflow. 

454 (5) Our results seem to confirm the conclusion of Moebius et al. [1987] and many other 

455 studies afterwards, in that acceleration processes in the near-Earth magnetosphere are 

456 mass dependent, because it is more effective for the heavier ions (0+/H+ ratio higher). 
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Appendix A: Construction of RAPID 0+/H+ ratio 

457 We shall refer to the corresponding differential particle fluxes, jj, from the eight en- 

458 ergy channels, {E^^ , E^~^), of the RAPID instrument as j^^ to j^^ and jf^ to j^^ , 

459 respectively (see Table 1). 

460 The lowest energy channel for oxygen, E^^ (with energy range 84-274 keV), is con- 

461 taminated by protons, and cannot be used [Daly and Kronberg, 2010]. The count rates 

462 in energy channels above 1 MeV (i.e., Eq, E-j and E^) are usually extremely low for both 

463 species. Useful O"*" channels are therefore energy channels E^^ - Ef^ , covering energies 

464 from 274 to ~955 keV. For protons, corresponding to this energy range, parts of E^'^ and 

465 E^^ have to be used. We have decided to use both these H+ channels rather than only 

466 one E^'^ and the part of the , respectively, for better statistics. To construct O"^ and 

467 H+ energy bins with matching energies, we proceed as follows: 

1. Determine the integrated oxygen intensity, J^^ , at energies from 274 up to ~955 
keV. The integrated intensity is determined as the differential fluxes, jf^ multiplied by 
respective A.E^'^: 

J^' = Ur ■ AEr + jT ■ + jT ■ + Jt ■ ^Ermicm^ ■ sr . .)]; 

2. Obtain spectral index 7, in order to cut the part of the E^'^ , see Figure 7 for 
illustration. For this we use fluxes from two adjacent energy channels {E^'^ and E^'^) 
and their effective energies, E^ff, the energy at which the spectrum, A ■ E~^j:, with given 
A and 7 has the differential flux equal J/AE: 

7 = Hjf/jf)/HE5u/E5"f,), 
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468 where Egff in this case is calculated as the geometric mean energy of the respective 

469 energy channels. See more details on these calculations in the RAPID Cahbration Report 

470 [Kronberg and Daly, 2009], chapter B.l; 

3. Determine j^^ — i.e., the H+ differential flux at the virtual energy E^^ , 274-374 
keV (which matches the 0+ energy channel), with effective energy -E^j^. = -\/274 • 374 
keV and determined 7, see Figure 7: 

■H+ ■H+ -i-ln(E"^, /E^^fA). 

471 4. Determine integrated proton intensity, at energies from 274 up to 962 keV: 

472 J^^ — j^^ ■ A.E^'^ + j^^ ■ AE^'^ . With this procedure, we estabhsh a virtual proton 

473 energy channel with the same lower energy threshold as for the oxygen. The slightly 

474 different upper energy level (962 and 948 keV, respectively) does not play any role here; 

475 5. Deflne the intensity ratio: ^+ — 

In order to compare our results with other observations we need to use the energy 

density ratio. To establish the energy density ratio, we assume that the relative geometric 

factors are correct for both species: 

0+ _ 7rV2^ -y/E^f ■ JO^AE _ a.jO^ 

nV2^.^j.J^-AE~ J"'' ^ ^ 

476 where 0+ = y/rno ■ \[EeJf ■ J^^ and — ^Jmu ■ \jE^j ■ J^^ . 

ill This calculation is based on the assumption that the effective energy is equal to the 

478 geometric mean of the corresponding energy thresholds. However, in our case the width of 

479 the energetic channel is quite large and this will lead to the deviation of the energy density 

480 from the value calculated using the effective energy as the geometric mean. The way to 

481 calculate this deviation one can find in the the RAPID Calibration Report [Kronberg and 
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482 Daly, 2009], chapter B.3. The deviation is estimated to be ~65% and error bar ±30% 

483 from the value of the energy density calculated in Equation Al. For these calculations 

484 the typical range of 7 values derived from our data base were taken: for O"*" 7 =2-3.5 and 

485 for H+ 7 =3.5-6.5. The statistical errors of the energy density and the error due to the 

486 large width of the energy channels are added in this case. 
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Figure 1. Maps of the location of the energetic oxygen observations, from left in 



XYgse, YZgse dayside {Xgse > 0) and YZgse nightside {Xgse < 0) planes. 
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Table 1. Lower thresholds of energy channels for protons and oxygen for RAPID. 

Note that energy channels for H+ and 0+ are not overlapping. To make the 0'''/H+ ratio 

meaningful, we therefore fit a spectrum and derive an additional H"*" energy channel for 

RAPID, see Appendix. 



Channel 


Energy range 


, keV Energy range 0+, keV 


1 


28-64^ 


84^ 


2 


75 


274 


3 


92 


414 


4 


160 


498 


5 


374 


638 


6 


962 


948 


7 


1885 


1414 


8 


c 


2539 


Upper 


4007 


4046 



^ Energy gap between hydrogen 1 & 2 
^ CNO channel 1 contaminated at times, suppressed 
8th channel not accessible for hydrogen 
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Table 2. Response of ion intensity to changes in tlie Dst index. 



Particles 


Compression^ 
Dst>15 nT^ 


Growth phase^ 
-30<Dst<15 nT 


Weak-moderate'^ 
-105<Dst<-30 nT 


Quiet-Storm time"^ 
-105<Dst<0 nT 


Min-Max*^ 


H+ ~10 keV 


no change 


no change 


+1.6x5 


+1.4x 


+2x 


0+ ~10 keV 


-2x8 


+1.5x 


+8x 


+14x 


+13x 


0+/H+ ~10 keV 


-2x 


+1.7x 


+5x 


+10x 


-F22x 


H+ >274 keV 


+2x 


+1.5x 


+1.6x 


+2.3x 


+2.4x 


0+ >274 keV 


-8x 


+3x 


+1.6x 


+3.6x 


+30x 


0+/H+ >274 keV 


-16x 


+2x 


no change 


+1.6x 


+30x 



^ Difference between ion intensities at ~22.5 nT and ~0 nT (here and further, the mean 

values of the corresponding bins), namely during magnetospheric compression. 
^ Difference between ion intensities at ~0 nT and ~-22.5 nT, namely between beginning 

and end of the growth phase. 
^ Difference between ion intensities at ~-37.5 nT and ~-97.5 nT, namely between weak 

and moderate storms. 

Difference between ion intensities at ~0 nT and ~-97.5 nT, namely between quiet and 

storm time. 

Difference between minimum and maximum ion intensities in Figure 2. 
^ This value and corresponding values in the next columns are the thresholds of the 

ranges. 

g a_*^)) -g (decrease in * times and "-|-*x" is increase in * times. 
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Dst, nT 

Figure 2. a) Dependencies on Dst index of (a) proton and oxygen intensities at ~10 
keV, (b) proton and oxygen intensities for >274 keV, (c) ratio O^/H"*" for intensities at 
~10 keV and (d) ratio 0+/H+ for intensities at >274 keV. The dashed hne in panel 3 
indicates a hnear fit to the ratio values. Error bars indicate 99% confidence intervals. 
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Figure 3. Dependencies on AE index of (a) proton and oxygen intensities at ~10 keV, 
(b) proton and oxygen intensities at >274 keV, (c) ratio 0'^/H+ for intensities at ~10 
keV and (d) ratio 0"''/H+ for intensities at >274 keV. The dashed line in panel 3 indicates 
a hnear fit to the ratio values. Error bars indicate 99% confidence intervals. 
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Figure 4. Dependencies on Clock Angle of (a) proton and oxygen intensities at ~10 
keV, (b) proton and oxygen intensities at >274 keV, (c) ratio 0+/H+ for intensities at 
~10 keV and (d) ratio 0+/H+ for intensities at >274 keV. 
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Figure 5. Dependencies on solar wind density of (a) proton and oxygen intensities at 
~10 keV, (b) proton and oxygen intensities at >274 keV, (c) ratio 0+/H+ for intensities 
at ~10 keV and (d) ratio 0+/H+ for intensities at >274 keV. 
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Figure 6. Dependencies on solar wind pressure of (a) proton and oxygen intensities at 

~10 keV, (b) proton and oxygen intensities at >274 keV, (c) ratio 0'^/H+ for intensities 
at ~10 keV and (d) ratio 0+/H+ for intensities at >274 keV. 
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Table 3. 0"^/H^ ratios of energy density depending on the disturbance level and the 
location. 



Satellite, 


Energy range, 


0'''/H"'", Energy density 




instrument 


keV 


Quiet time 


Disturbed time 


Cluster/RAPID 


274 keV - ~955 keV 


0.22±0.12% 0.72±0.44'^ 


1.04±1.33% 1.11±1.35^ 


Cluster/CIS 


~10 keV 


0.038±0.0081% 0.034±0.0031^ 


0.083±0.038^, 0.2±0.12^ 


AMPTE/CCE'= 


1-310 keV 


0.03 


0.34 


AMPTE/CCE'i 


1-310 keV 


0.01 


0.61 


Geotail/EPIC^ 


9-210 keV 


0.05-0.1 


0.2-0.6 



^ Quiet (AE<100 nT) and disturbed conditions (AE>300 nT) are based on the AE 

index. 

^ Quiet (Dst~ nT) and disturbed (Dst is between -100 and -30 nT) are determined 

based on the Dst index. 
^ Measurements taken from Gloeckler et al. [1985]. 
^ Measurements taken from Hamilton et al. [1988]. 

^ Measurements taken from Nose et al. [2001]. 
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Figure 7. Bringing RAPID H+ channels to the same energy range as O"^ channels: 
the differential flux in channels E^^ and E^^ is fltted to a power law, that flt is used to 
determine the differential flux, j^, in the virtual channel E^. E^ff is the energy at which 
the spectrum with given A and 7 has the differential flux equal J/AE, where J is the 
integrated intensity. 
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